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A number of Galactic sources emit GeV-TeV gamma rays that are produced through leptonic 
and/or hadronic mechanisms. Spectral analysis in this energy range is crucial in order to un¬ 
derstand the emission mechanisms. The HAWC Gamma-Ray Observatory, with a large field of 
view and location at 19° N latitude, is surveying the Galactic Plane from high Galactic longitudes 
down to near the Galactic Center. Data taken with partially-constructed HAWC array in 2013- 
2014 exhibit TeV gamma-ray emission along the Galactic Plane. A high-level analysis likelihood 
framework for HAWC, also presented at this meeting, has been developed concurrently with the 
Multi-Mission Maximum Likelihood (3ML) architecture to deconvolve the Galactic sources and 
to perform multi-instrument analysis. It has been tested on early HAWC data and the same method 
will be applied on HAWC data with the full array. I will present preliminary results on Galactic 
sources from TeV observations with HAWC and from joint analysis on Fermi and HAWC data in 
GeV-TeV energy range. 
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1. Introduction 

GeV-TeV gamma rays are produced by energetic particles interacting with other particles or 
magnetic field. Supernova remnants (SNRs) and pulsar wind nebulae (PWNe) are most common 
Galactic sources in the TeV energy range. In PWNe, electrons are accelerated by the termination 
shock where the pulsar wind collides with the surrounding gas, emitting gamma rays via inverse 
Compton scattering. In SNRs, which are postulated as the cosmic-ray acceleration site, hadrons 
gain energy via diffuse shock acceleration, emitting gamma rays via pion decay. About half of TeV 
Galactic sources are still unidentified (UID) wifhouf associafing in ofher wavelengfhs. Observafions 
of GeV-TeV gamma rays are the key to distinguish between leptonic and hadronic processes which 
will aid in source identification. 

The inner Galaxy region contains the strongest TeV Galactic sources other than the Crab Neb¬ 
ula. H.E.S.S. imaging atmospheric Cherenkov telescopes (lACTs) [1, 2] surveyed the Galactic 
plane in the southern sky and air shower arrays Milagro [3] and ARGO [4] surveyed the Galactic 
plane in the northern sky. Targeted observations have been performed by lACTs H.E.S.S., VERI¬ 
TAS, and MAGIC (for example, [5, 6]). Fermi Earge Area Telescope (Eenni-EAT) performed an 
unbiased all sky survey in the energy range of 20 MeV-300 GeV. Their third source catalog [7] is 
published based on the first four years of science data, and another catalog of sources above 10 GeV 
is published with the first three years of Fermi-LP^I data [8]. 

The High Altitude Water Cherenkov (HAWC) Gamma-Ray Observatory, located at Sierra 
Negra, Mexico at 4100 m a.s.L, is sensitive to gamma rays in the energy range of 100 GeV to 
100TeV [9]. HAWC consists of an array of 300 water Cherenkov detectors (WCDs), covering an 
area of 22,000 m^. Each WCD is 7.3 m in diameter, 4.5 m in depth, and filled with ~ 200,000 E of 
purified wafer. Eour upward-facing phofomulfiplier fubes (PMTs) are affached fo fhe boffom of each 
WCD fo defecf the Cherenkov radiation produced by the secondary particles in an air shower. One 
high-quantum efficiency 10-inch Hamamafsu R7081-MOD PMT is in fhe cenfer of each WCD and 
fhree 8-inch Hamamafsu R5912 PMTs are af 1.8 m from the center [10]. The HAWC observatory 
has a high duty cycle of > 95% and a large instantaneous field of view of 2 sr. 

A major challenge analyzing the emission from the inner Galaxy region, where has high po¬ 
tential of source confusion, is to deconvolute and identify sources due to limited angular resolution 
of TeV air shower arrays. A likelihood framework [11] has been developed to simultaneously fit the 
positions and fluxes of mulfiple sources, allowing fo defermine fhe number of sources in a region 
of inferesf (ROI). In fhese proceedings we presenf a TeV gamma-ray survey of fhe inner Galaxy 
region of Z G [-1-15°, -|-50°] and b G [-4°, -|-4°] using fhe dafa faken with a partially-completed 
HAWC array, followed by a joint analysis using Fermi-LAJT data. 

2. Data and Analysis 

The full HAWC array was completed in March 2015, but science operations already started in 
August 2013 with a partially-built array. In these proceedings, the analysis is performed on data 
taken with the partially-constructed HAWC array, hereafter referred to as HAWC-111, between 
August 2, 2013 and July 8, 2014. During this construction phase, the active WCDs grew from 
108 (362 PMTs) to 134 (491 PMTs). Data taking was occasionally interrupted during this period 
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because the deployment of the detectors took precedence. The dataset contains 275 ± 1 source 
transits in the inner Galaxy region. 

The energy of a primary gamma-ray particle is correlated to the size of the air shower. Data 
are binned in 10 bins according to the fraction/ of PMTs that are triggered in an air shower with 
respect to the total live PMTs. HAWC-111 detects air shower events at a rate of 15 kHz. However, 
> 99% of air showers that trigger the detector are induced by charged cosmic rays, which are the 
main background for gamma-ray observations. A series of cuts are applied to the data [12] to 
remove cosmic ray events. Cuts are optimized in each / bin to maximize the significance of the 
Crab nebula, which is the brightest steady source and a standard candle in the TeV energy range. 

The arrival direction of each air shower event is reconstructed using the relative time of the 
PMT signals when the air shower crosses the array [13]. The arrival directions are then filled info 
a signal map using HEALPIX pixelafion scheme [14]. The average pixel scale is 0.11°, which is 
smaller fhan fhe angular resolufion of fhe HAWC- 111 defecfor. The background is esfimafed from 
fhe dafa using fhe mefhod of direcf integration [15]. The infegrafion lime of fwo hours is chosen lo 
emphasize fhe sources smaller fhan 30° in angular size. 

A sel of Monte Carlo simulations, which simulale fhe defector response lo gamma rays and 
cosmic rays, are used in fhe maximum likelihood mefhod. These simulations are generated using 
fhe CORSIKA package [16], which simulates fhe shower propagalion in fhe almosphere and fhe 
GEANT4 package [17], which simulates fhe response of WCDs to air shower particles. 

2.1 The Likelihood Method 

The likelihood mefhod begins wilh building a source model. A poinl source is characlerized 
by ils position and speclrum. The position is described by righl ascension and declinalion. The 
speclrum is given by a simple power law 



( 2 . 1 ) 


where Iq is fhe flux normalization, Eq is fhe pivof energy, and a is fhe specfral index. The sludy 


presenfed here assumes a fixed index of 2.3 according lo fhe average measured specfral index for 


known Galactic objecls [1]. 

In each pixel of each/ bin, fhe observed number of evenls follows a Poisson disfribulion. The 


mean of Ihis Poisson disfribulion is fhe expecfed signal A, which is fhe number of background 


evenfs b in fhe dafa plus fhe expecfed number of gamma-ray evenls y given a source model con¬ 
volved wilh fhe defecfor response from fhe detector simulalions. The probabilify of observing a 
signal N given an expecfed signal A is 



( 2 . 2 ) 


The log likelihood of having fhe observalion given a paramefer sef 0 in fhe source model is 
fhe sum of fhe log likelihood of each pixel in an ROI and in each/ bin. 


fbinsROI 


iog^{e\N)= £ ^{NijiogXij-iij) 


(2.3) 


i j 
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where Nij and Xij are the observed and expected signal in the 7 th pixel of the ithf bin, respectively. 
The term Nij ! is discarded since it is given in the observation and independent of the parameters in 
the source model. The MINUIT package [18] is used to maximize the log likelihood by varying the 
parameter set 6 in the source model, i.e. to find the most likely source model given the observation. 

The size of ROI is chosen to be larger than the angular resolution of the detector and possible 
source extent to include more photons from a given source. However, there is a high potential of 
source confusion for an analysis of the Galactic plane in the HAWC-111 data. In other words, it 
is not always possible to choose an ROI that contains only one source without contamination by 
photons from another source. In this case the source model may need to build with more than one 
source and the expected signal becomes 


j ~ y ^ > (2-4) 

k 

where yijk is the expected number of gamma ray events from the kth source in the jth pixel of ithf 
bin and bij is the background in the jth pixel of ithf bin. 

To compare different source models, a likelihood ratio test is performed. Given the log likeli¬ 
hood of the background-only model log^o and the one-source model log^\. A test statistic {TS) 
as defined in 


TS=-2{log^o-log^i) (2.5) 

indicates how much the one-source model is better over the background-only model. TS follows a 
distribution with the degrees of freedom (DoF) being equal to the difference on the number of 
free parameters between two models according to Wilks’ theorem. Then TS can be converted to a 
p-value of at what level the observation is consistent with the background-only hypothesis. 

The significance map of the inner Galaxy region as shown in Fig. 1 is made by moving a pu¬ 
tative point source through each pixel, performing a maximum likelihood fit on flux normalization 
with spectral index fixed at 2.3, and converting the TS value to significance according to Wilks’ 
theorem. 

Due to limited sensitivity and angular resolution in this data set, all sources in the model are 
considered as point sources. To reduce the number of sources and number of free parameters in the 
likelihood fit, five ROIs are chosen in the inner Galaxy region covering all the pixels with > 5a 
pre-trials significance. The position and size of each ROI, as shown in Fig. 1, is carefully selected 
to avoid a source at the boundary. To decide the number of independent sources that are needed 
to properly model each ROI, an iterative likelihood ratio test is performed, in which an additional 
source characterized by three free parameters, right ascension, declination, and flux normalization, 
is added into the model at each step. This iterative procedure is repeated until a priori chosen 
threshold of NTS < 15 (~ 3 a with 3 DoF) with an additional source. This procedure provides the 
number of sources in each ROI, and results the position and an initial estimation of the flux of each 
source candidate. After the iterative procedure is done for each ROI, a global fit on source fluxes 
in all regions is performed with fixed source positions to minimize photon contamination from 
sources just outside of a ROI. Finally, a single-source likelihood fit is performed on each source 
candidate by treating other source candidates as part of the background to obtain the flux and TS. 
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significance 


Figure 1: Significance map of the inner Galaxy region for 283 days of HAWC-111 data. The open squares 
mark TeV sources in the TeVCat [19] and the open circles indicate GeV sources in Fermi-hPJF 3FGL [7]. 
Five black circles show the five ROIs used in this analysis. 


3. Results 

Ten source candidates are identified using the likelihood analysis with > 3a post-trials, with 
eight of them with tentative TeV associations. More discussion on the TeV counterparts can be 
found in [20]. As shown in Fig. 1, a number of these source candidates are spatially coincident 
with GeV sources in Ferm/-LAT 3FGL or Fenni-LKF IFHL catalog. However, due to the lim¬ 
itation of angular resolution in HAWC-111 data, it is not clear if the TeV and GeV sources are 
associated without further investigation. The results from a joint likelihood analysis using 3ML 
[21] framework on source spectra in GeV-TeV energy range with HAWC-111 and Fermi-LPJF data 
will be presented at the ICRC. 
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